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centre of gravity ft/sec 
velocity 


wind velocity ft/sec 
initial wind speed ft/sec 
range £E 
crossrange ft 
height fats 
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I. INTRODUCTION 

A. ARMT OPERATIONAL CONCEPT 

Conventional artillery uses unguided projectiles to 
destroy a variety of targets including armored vehicles. 
Ballistic information and professional Sse Sire used 
to direct fire onto a target with little chance of a4 direce 
hit in the first few rounds. A forward observer provides 
feedback to the guns for correction until the target is 
destroyed with subsequent rounds.Generally this approach is 
limited to the 30 kilometer range and often point target 
locations are not Known. 

The ARMT projectile,formerly called AIFS, changes the 
conventional concept by providing a fire and forget,boosted 
and terminally guided projectile with a highly predictable 
trajectory.Conventional heat seekers and target identifica- 
tion systems are short range and require a blunt nose which 
has high drag. Active radar is not practical. Long ranges of 
60 kilometers are achieved by using solid fuel ramjet 
propulsion with an axisymmetric air inlet spike. An 
accurately predicted trajectory ig achieved by ensuring 
thrust equals drag to produce zero axial acceleration and 
create a pseudo-vacuum trajectory during powered flight. At 
a preset time (near ideal condition ramjet burnout) the 


inlet sapike is discarded and canards are deployed for 
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terminal guidance control. A passive or reflected target 
energy system will be used for terminal guidance. Research 
is ongoing and sponsored by DARPA. Projectile 


characteristics are at Appendix A. 
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Figurerdi-i-aARMT Trajectory Concept 


B. THRUST DRAG BALANCE CONCEPT 

Balanced thrust drag implies that the thrust equals the 
drag vectors during powered flight. The aim of the balance 
is to negate meteorological variations throughout the 
trajectory by reducing all force vectors to zero except 
gravity. The resulting pseudo-vacuum trajectory is highly 
predictable. The balance is accomplished by measuring axial 
acceleration and maintaining it at zero.Appendix B explores 
this concept in detail. The potential errors include 
transients due to wind induced reactions and measurement 
errors in the technical application of the balance in 
flight. Further, any asymmetries such as aerodynamic 
misalignment’ or thrust will produce a lateral force which 
must be offset by a smali roll rate(sixth degree of 
freedom). This allows longitudinal forces to be the 
significant forces. 

The concept of thrust drag balance is best demonstrated 
in Figures 1-2 and 1-3.In Figure 1-2 the centre of gravity 
of the projectile is flying at an angle of attack (alpha) 
with a velocity vector (Vcg) relative to a fixed earth 
reference frame. The aerodynamic forces are lift(L) and 
drag(D). Assuming static static stability and neglecting 


transient dynamics, the angle of attack becomes zero, which 
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Reauces, |) Litt to zero.Tf thrust equals drag then only the 
gravitational force exists. 

If a crosswind(Vw) is applied, the projectile cocks 
into the wind because of the aerodynamic force acting on the 
centre of pressure. Adjusting the thrust and again assuming 
static stability the ena can be offset as shown in Figure 
1-3 by applying the correct thrust.For future reference the 


component of drag normal to the plane containing Vcg is Dw. 
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Figure 1-2 Projectile Dynamics Thrust Equals Drag 
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Figure 1-3 Projectile Dynamics With Crosswind 


If a tailwind or headwind is applied, thrust is adjusted 
to compensate in the same manner. 

The dynamics or transient effects due to wind which 
create damped oscillatory motion have a cumulative net lift 
effect( as shown in Figure 1-4) as angle of attack is 


proportional to lift.The area under the curve in Figure 1-4 


does not average to zero due to the damping. For future 


reference this lift due to wind will be referred to asa Lw. 


Figure i-4 Angle Of Attack Versus Time 


Calculations integrating alpha(t) twice to demonstrate the 


above in the yaw or pitch plane are at Appendix B. 


See bHeE StS SCOPE 

It is emphasized that the scope of this thesis is to 
analyze dispersion from the pseudo-vacuum trajectory at 
terminal guidance time for the eight inch ARMT projectile.In 


colloquial terms the “basket from whveh): the projectile 
will be guided to target was determined. Further, this 


basket is based on thrust errors and wind effects for a 


LS 


launch angle of 45 degrees and muzzle velocity of MACH 2.2 
at sea level. Linearized aerodynamics were used in all 


calculations. No attempt was made to establish a "footprint" 


for projectile impact. 
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II. TRAJECTORY ANALYSIS METHODOLOGY 

A. VACUUM TRAJECTORY ANALYSIS 

Deo snort ~and” simple" BASIC “proegram was written to 
establish the vacuum trajectory and fuel burn time of the 
projectile in a two dimensional scenario where thrust 
exactly equaled drag. The Standard U.S. Atmosphere was used 
and the drag coefficient linearized {Reference 1]. The 
program listing and baseline data are at Appendix D. The key 
parameters and equations used are listed at Appendices A and 
D.The baseline pseudo-vacuum trajectory is shown in Figure 


Z2-1. 


B. THREE-DEGREES-OF-FREEDOM MODEL 

The computer simulation was established using the 
Interactive Ordinary Die sonentwas Equation (IODE) library 
program at NPGS. A listing of the three degrees of freedom 
equations is in the IODE specifications at Appendix E. The 
Simulation is based on wind tunnel results from Norden 
Systems (References 2 and 3] and assumes a horizontal plane 
at 30000 feet altitude.This height was chosen as significant 
frrgnt time is at altitude and the time of flight in this 
plane was near actual flight time. Primarily the model was 
used to study crosswind effects but resonance with the 


projectile short period was simulated. Different thrust 
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scenarios 


the X-Y plane. 


*10° 


HEIGHT (FEET) 


0.0: 72.0¢% 4.0; 26.0 


Figure 2-1 


were also used. Figure 2-2 depicts the model 


burnout time=101 secs 
burnout range=175100 feet 
burnout height=10990 feet 
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ARMT Pseudo-vacuum Trajectory 
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Figure 2-2 Three Degrees Freedom Model (Yaw Plane) 


The crosswinds were varied in intensity and duration to 
analyze crossrange dispersion. As well thrust was varied as 
a function of drag by scaling the drag. Cmq predictions were 
approximately -190 to -205 by Norden Systems [Reference 2] 
and an analysis was done to determine the effects of varying 
this damping coefficient. To further study the projectile 
characteristics, trajectories were simulated with no thrust 


fon Wass aurenuation with varying Wind profites: Figure 2-3, 
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shows the baseline three-degrees-of-freedom trajectory 


thrust equals drag with no wind. 
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Figure 2-3 Baseline Three-Degrees-Of-Freedom Trajectory 
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C. FIVE-DEGREES-OF-FREEDOM MODEL 

Again the IODE program was used to solve a five degrees 
of freedom system of equations. The roll plane was neglected 
for the following reasons: 

-the program could not accept sufficient met aena to 
define six degrees of freedom. 

-the fins are expected to damp out initial spin due to 
guns rafling ©9000, RPM). in. the first few seconds of flight. 
Sufficient spin to average projectile asymmetry errors to 
zero but not enough to significantly affect dynamics will 
remain.(ie. gyroscopic effects and Magnus effects) 

-the initial spin damping dynamics are beyond the scope 

of this thesis. 
The projectile was assumed to be axisymmetric allowing for 
ease of calculations in the yaw and pitch planes. Figures 
a2 and 2-4 depict “the yaw “and pitch planes for the 
model,respectively. 

a Listing of the” “equations ss provided ‘in the IGDE 
specifications at Appendix F. Trajectory simulations were 
run with varying wind profiles and thrust error scenarios. 
Terminal guidance time (Tc) was determined by the baseline 
five degrees of freedom trajectory shown at Figure 2-1 as 
101 seconds. This trajectory has no wind and thrust exactly 


edualted” drag “until” disceard of the air "inlet " spike and 
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eanard deployment.All simulations assumed a 45 degree launch 
angle and Mach 2.2 muzzle velocity. Dispersion from the 
pseudo-vacuum trajectory was measured at terminal guidance 
time including the following variables: 

-X=range at 101 seconds 

-Y=crossrange at 101 seconds 

-Z=height at 101 seconds 

-Tb=burnout time 

-Vcg=speed at time 101 seconds 


Dispersion is defined as the difference between value for 


perturbed trajectories and the baseline trajectory. 
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Figure 2-4 Five-Degrees-Of-Freedom Pitch Plane Model 
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Di. SIX-DEGREES-OF-FREEDOM MODEL 

The Ballistics Research Laboratories (Aberdeen) six 
degrees of freedom model HTRAJ was obtained and adapted for 
use on the NPGS IBM 370 system. The adapted listing is 
included as Appendix G. Program compilation was achieved but 
output data was not satisfactory due primarily to input data 
syntax problems. Further,the program requires extensive 
modification to be able to vary thrust as a function of drag 
with the ease of the previous models already described. Time 
was not available to correct these problems. Detailed 


documentation on the program HTRAJ is available through BRL. 
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Ifl. TRAJECTORY RESULTS 


Figures 3-1 ,and 3-2.,show the effects of . constant 
thrust errors as a percentage of drag with no wind.Recall 
that the three-degrees-of freedom trajectory is at constant 
altitude.The figures therefore represent a thrust error 
sensitivity in terms of burn time and range variations due 
to drag or thrust.The trajectories all pass nearly through 
the 218000 foot range mark at t=100 seconds, but the burnout 
ranges are significantly different due to the varying burn 
times resulting from projectile thrust.Figure s-2 
demonstrates the obvious conclusion that the burnout time 
decreases as thrust increases. 

Figure 3-3 combines the effects of constant 
crosswinds and thrust errors to demonstrate resultant cross 
range dispersion.The constant wind velocity lines are very 
nearly linear, and the pattern is essentially a shift about 
the zero-zero point. If thrust is leas than drag, the 
dispersion ia in the direction of the wind vector.If thrust 
ig greater than drag then the overcompensation produces a 


dispersion opposite to the wind vector. 
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Figure 3-1 Thrust Error Sensitivity (Range) 
Three Degrees Of Freedom 
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Figure 3-2 Thrust Error Sensitivity (Burntime) 
Three Degrees Of Freedom 
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Figure 3-3 Crosswind Sensitivity Three Degrees Of 


Freedom 


Of interest was the dispersion due to the same wind 
profiles without thrust or mass loss. This represents an 
absolute dispersion with a given wind profile. These results 


are in Table 3.1. 


Table 3.1 No Thrust Crosswind Effects 
Three Degrees of Freedom 


Thrust=O t=125 seconds 


vw (ft/sec) dispersion (feet) 
crossrange 
‘e) fe) 
60 L935\.76 
30 293.80 
120 392.02 


The data clearly shows that crosswinds have a 
minimal effect in the thrusted scenarios.The dispersion in 
the thrust-equals- drag case is due to the transients (Lw) 


and is negligible as well. This is illustrated at Table 3.2. 
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Table 3.2 Dispersion Due to Net Lift 
from Transients 


Thrust=drag 


vw (feet/sec) dispersion (feet) 
crossrange 
re) fe) 
60 sons 
90 O27 
120 -036 
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Angle of attack damped to virtually zero during the 
trajectories (see Figure E-3) but the absolute maximun 
values were also recorded.Figure 3-4 demonstrates the 
relation between angle of attack (alpha) and the wind 
velocity.Thrust errors had virtually no effect on the angle 


of attack magnitude. 
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Figure 3-4 Wind Effects On Maximum Angle Of Attack 


2. Damping Coefficient Effects 
The magnitude of Cmq was determined to be 
approximately -200 in the flight speed ranges of the ARMT 
projectile (Reference 4]. To determine if this coefficient 
Significantly affected dispersion Cmq was varied from O to 
1.9 x 106 and the dispersion recorded. Figure 3-5 shows the 


relation between Cmq and dispersion (crossrange) with a 


constant wind of 60 feet/second and time of 125 seconds. 
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Figure 3-5 Damping Coefficient Effects 
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Careful examination of Figure 3-5 shows that dispe- 
rsion increases Slightly for realistic values of Cmq and is 
therefore proportional to Cmq.Due to the small sensitivity 
of dispersion to Cmq for all practical purposes one can 
conclude dispersion is independent to Cmq.(The increased 
dispersion proportional to Cmq can be explained as the 
difference in the areas under a highly damped sinusoid 


versus a non-damped sinusoid is considerable.) 


3. Sinusoidal Wind Analysis 

To test the effects of phase angles and wave- 
lengths on the projectile dynamics a sinusoidal wind was 
input to test the model.The wind (Vw) was of the form Vo 
Sint2 PI Vo t/lambda + phase) where Vo was set at 60 
feet/second.Figures 3-6 and 3-7 show the trajectories for 
the sine and cosine winds to burnout at thrust equal to 90% 
drag.(Lambda=range and lambda=.5 range)Of note is the fact 
that the sine wind produces a continuous dispersion while 
the cosine wind does not when thrust does not equal 
drag.This discovery led to an expanded study of different 
phase and wavelength sinusoidal wind inputs for various 
thrust scenarios.As well a simplified model was constructed 
as shown at Figure 3-8 and calculations confirmed the three 


degrees of freedom model results.Figures 3-9 and 3-10 depict 


the results of the calculations. 


QJ 
m 


DISPERSION (FEET) 


8.0 


Ue? 


6.4 


LEGEND 
LAMBDA=304000 


5.6 


4.8 


4.0 


See 


2.4 


1.6 


0.8 


0.0 


0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 
TIME (SECONDS) 


Figure 3-6 Sinusoidal Wind Effects (Sine) 
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Figure 3-7 Sinusoidal Wind Effects (Cosine) 
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Figure 3-8 Simplified Three Degrees Of Freedom Model 


The Simplified model was used to analytically 
evaluate the crossrange dispersion (y) in the cases of sine 
and cosine winds to confirm the validity of the three- 
degrees-of-freedom model.Calculations follow: 

General equations of motion in the y direction from 


Figure 3-8 are: 
M dVy/dt =Dy which approximately equals D(Vw-Vy)/Vxo 
if Vxo=Vo=constant.Then transforming the equation we get: 


dVy/dt + D/(MVo)Vy = D/C(MVoO) Vw lt) 


ie 
w 


or in time differentiated notation: 
y+ By = B Vuct>) is D/Civom—e 


If a cosine wind Vw(t)=Vwo cos (2 PI Vo t/A) is 


and the differential equation solved we have: 


y(t)/% =K }e-Bt -cos(2 PI Vo t/%) + 


DAsin(2 PI Vo t/A) 


2 PI m Vo2 


let K= 2 $172 me@Vozl Vwo 
DA Vo 

1 +4 PT DA 2 

Ge m || 


applied 


Figure 3-9 graphically depicts the results as t increases.A 


zero asymptote is observed and this conforms with 


results. 


model 


Now if a sine wind Vw(t)= Vwo sin (2 Pi Vo t/A) is applied 


and the differential equation solved: 


y(to/A =K Div (l-eos(2 PI Vo E72.) 
4 PI 1/2 m Vo2 


sin¢2 PI VO t/2X) 44 PECI/2 mVoZ) C1—e7-BeD 


DA 


Figure 3-10. graphically « depicts», these “results 
increases.The asymptote is 


K. | (4. PT 172 moVoZzZ) Ci-esBbe) 
D> 


and this conforms to the model also. 
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Figure 3-9 Cosine Wind Analysis Results 
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Figure 3-10 Sine Wind Analysis Results 
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The constant wind case was also derived: 


Vwtt) Vw andthe solved differential 
equation is: 


VGED=EVekt #=hVG maVvo7 Die > Bt 


and as t increases y approaches Vw t. 


Specific calculations were done to confirm 
dispersion results at burnout.The calculations were based 


on: 


y=dy/dx X = 1/Vo dy/dt X 

Vwo=60 ft/sec 

Vo=2190 ft/sec 

~ =304000 ft 

t= A~/Vo 

D=1/2 9 Vo2 A Cd (=226 lbf) 
Thrust=.9D therefore D=22.6 lbf 


M=2235 slbm 


Taking the time derivative of the cosine wind y(t) 


dispersion formula derived earlier we get 


dy/dt=K|-Be-Bt + 2 PI Vo sin(2 PI Vo t/a) + 


A 
DA Mec eoneces tz PE Vo E/) 
4 PI 1/2 mVo2 oN 


at t=A/sVo the equation becomes: 
dy/dt = K(-Be-BA/Vo + B) 


and substituting the values into the approximationn 
for ys) y=O) BE. 
Now, performing the same operations on the sine 


wind dispersion equation 


dy/dt=K DA , §2 PT) Vosesiné2) Pigvese/a) 
4 PI 1/2 mVo2 oN 


-2 PI Vo cos(2_ PI Vo t) +(4 PI 1/2 m Vo), 
mS x DA 


(Be~Bt) 
and at t=WVo 
dy/dt, =K:| -2. PI Vo, -4. PR 1/2 mavo2sbaslE 
» DX» 
and y=20 ft approximately. 


Final iy; taking the constant wind dispersion 


formula and taking its derivative 
dy/dt = Vw- Vw e Bt 


and y= 10 ft approximately. 
These results are consistent with the three-degrees-of- 
freedom model considering the initial approximations. 
To study the short period resonance phenomenom the 


Phase and wavelength wind variations were documented in 


ageteait for a condition of thrust equal to 90% drag.Figures 
3-11 to 3-13 show the crossrange dispersion effects from 
wavelengths of S000 ¢( lambda = JORG range) ES 
303000(lambda=range) feet and the phase angle variations 
from O to 2 PI.At approximately 9850 feet the dispersion 
characteristics were no longer symmetric but remained under 
one foot.This is demonstrated at Figure 3-1¢.A crossrange 
dispersion contour was constructed at Figure 3-15 but this 
did not emphasize the resonance phenomenon sufficiently due 
to linearized aerodynamics and small dispersion 
magnitudes.The parameters that varied most significantly 
near resonance were angle of attack and burn time.Indeed,it 
is clear that the projectile is unstable in this region with 
a maximum angle of attack of 1.21 radians(70 degrees) .Due to 
the increased angle of attack, drag increased and burn time 
decreased in the order of seven seconds (5%) compared to the 
baseline trajectory. This clearly affected the range at 
burnout .Due to large angle of attack the linearized 
aerodynamics are not valid.Projectile performance at wind 
resonance is worthy of additional analysis.It was decided to 
demonstrate resonance by plotting the maximum angle of 
attack (alpha) versus wavelength with different thrust/ drag 
scenarios.The resonance results are at Figure 3-16.The 
projectile natural frequency is estimated at 0.22 Hz. 


Resonance was confirmed in the 0.22 Hz regime for thrust 


errors under 10% including varying Vwo.This was necessary as 
the system is non-linear and Vwo could have affected 
resonance.It was noted that the curve shifted and increased 
or decreased in magnitude slightly in each case.This is due 
to the non-linearity of the system equations and the 
variations in Vo during flight.The absolute case ef) Jna 
thrust showed resonance at 10000 feet with a maximum angle 


of attack of 1.25 radians(72 degrees). 
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Figure 3-11 Dispersion Versus Wind Phase Angle(Large Lambda) 
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Figure 3-12 Dispersion Versus Phase Angle 


(Medium Lambda) 
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Figure 3-13 Dispersion Versus Phase Angle(Small Lambda) 
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Figure 3-14 Dispersion Versus Phase Angle 
Resonance Lambda 
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Figure 3-15 Phase Angle Versus Wavelength 
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Figure 3-16 Detrimental Resonance Between Projectile 
Short Period and Spatial Frequency of 


Wind Profile 


It is clear that if the wind profile experienced by 
the projectile has harmonic content near the natural 
frequency then the round will be unstable and/or have large 
dispersion at burnout. This is best illustrated by viewing 
the wind as Vw(t) with reference to the projectile. By a 
Fourier transform of the wind profile resonance with 
projectile natural frequency (fn) will be evident in the 
frequency domain.The wind oscillations need only be velocity 
variations with the appropriate harmonic content.Figure 3- 
17, which illustrates this concept, shows fn coinciding with 
@ maximum of Vw(f). One would anticipate stability problems 


for this case. 
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Figure 3-17 Fourier Transform Of Wind Profile Concept 
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Calculations for the long period resonance or 
phugoid show that the phugoid is irrelevant as flight time 
ia) too short: 

T=4.44*Vo0/g=301.8 seconds (Reference 3] 

The results of the sinusoidal winds were also 
compared with the aerodynamic jump phenomenom outlined in 
BRE-* report 1077 tReterence= 4° °)]" “At ’’ first “glance “the 
Sinusoidal wind results appear consistent with aerojump,but 
detailed calculations show otherwise. The rate of change of 
mora. “angle “of” attack’ is” critical’ to”“aerojump “and is 
easily incorporated into the aerojump equation.The results 
for the sine wind are a function of wavelength and of the 
same form as the simplified model results but the cosine 
wind and constant wind both yielded a small aerojump 
proportional “to the initial wind.The magnitudes of the 
aerojump dispersions were more consistent with the net lift 
transient effects. Further, continuous perturbations in the 
form of wind after time zero are not in the aerojump 
equations and for all practical purposes the phenomenom is 
not observed in the model.These results should have _ been 
apparent before calculations were performed. Aerojump 
equations are not appropriate for non-zero wind 
conditions.Clearly, when wind is included the aerojump 


equations cannot predict dispersion. 
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B... FIVE _DEGREES_OF FREEDOM 

As previously stated dispersion was measured from the 
ideal pseudo-vacuum trajectory at. a fixedy tame Lor 
initiation. of ..terminal guidancesstic) 101.» seconds. jhause 
error sensitivity in the zero wind condition is presented in 
Figures 3-18 and 3-19.In these two figures it is clear that 
thrust errors can have a significant impact on the location 
of ; the projectile. ar Phe lc higures so12 indicates that as 
Ehnusc increases the dispersion moves downrange as 
expected.Figure 3-19 shows that as thrust increases the 
height of the projectile at Tc decreases. Due to added 
thrust the projectile is located at a point further along 
the parabolic trajectory; hence the height is less. 

As well Figure 3-20 shows these effects on available 
burn time of the ramjet.For thrust equals drag burnout at 
101 seconds exactly equals the Tc.With less thrust, terminal 
guidance is initiated before all the propellant is burned. 
Conversely, when thruat exceeéeda drag, burnout oceursa peter 
to initiation of terminal guidance.Assuming this fixed time 
for initiation of terminal guidance, the projectile may 
experience up to seven seconds of unpowered flight (¢( 7% 
error) before the guidance ia deployed. The period of 
unpowered flight with associated high drag will cause major 


dispersions as shown in subsequent figures. 


OOSIZI OO0OO9LT OOSSLT O0OGZT OOST AT OOOFZLI 
(LAqd) FONVY 


5 10 
THRUST ERROR (%DRAG) 


0 


=o 


10 


| l 
OOSEZIT OOOEAT 


18 Thrust Error Effects (Range) 


Figure 3 


ip 
iB 


OSZTI OOcIT OSTTLT OOTTT OSOTT OOOTT OS6BOT OOG6OT OSssot OoogodT O0SZODOZOT 
(LAG) LHD 


+) 


THRUST ERROR (%DRAG) 


0 


=o 


10 


19 Thrust Error Effects (Height) 


Figure 3- 


BURN TIME (SECONDS) 


G8: 992 100 


97 


-10 —5 0 5 10 
THRUST ERROR (%DRAG) 


Figure 3-20 Thrust Error Effects (Burn Time) 


Crosswind sensitivity results were very much the same 
as the three degree sSimulation.In fact the results appear 
even more linear with the same rotation about the zero zero 
point as in Figure 3-21.The gross features of Figure 3-3 and 
3-21 are identical but Figure 3-3 shows dispersion at 
burnout whereas Figure 3-21 shows dispersion at Tc. 

The crossrange dispersion results of a no thrust 
projectile with constant crosswind profiles and with thrust 
equal to drag (‘transients effects) are at Tapes 3.3 and Sac 


respectively. 


Table 3.3 No Thrust Wind Effects 
Five Degrees 


mem we ewww we ae ae ae ae i es i i ea i es a es a es a a ee 


Thrust =0 time=101 seconds 
Vwxy (feet/sec) diapersion (feet) 
crossrange 
fe) fe) 
30 S21 
60 USSi5 
30 PLUS Se) ®) 
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Figure 3-21 Crosswind Effects Five Degrees 
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Table 3.4 Dispersion Due to Net Lift 
from Transients 


eee ae ae ae i a ei ae a ei ia iw a i a a we wa w= = i a a a a Se ee SS 


Thrust=drag time=101 secs 


Vwxy (feet/sec) dispersion(feet) 
fe) fe) 
30 »0202 
60 -0413 
930 | -0643 
120 20949 


Burn time was not significantly affected by crosswinds. 
Headwind and tailwind effects on range,height and 
burntime in different thrust scenarios are shown in Figures 
3-22 to 3-26. It is clear that dispersion from the) pseudo-— 
vacuum trajectory is minimized when thrust nearly equals 
drag.It is also apparent that the overthrusted trajectories 
involve greater dispersion as burnout occurs prior to 
terminal guidance and unpowered flight occurs for some time 
dependent on the scenario.This same error is observed in the 
case of headwinds and is most clear in the Figures 3-22B to 
3-25B about the thrust equal to drag axis. (Note that Figure 
3-22 to 3-25 inclusive have A,B and C sub-figures whieh are 
simply enlargements of key portions of the respective main 


figure). 
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Figure 3-22 Headwind Sensitivity 
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Figure 3-22A Headwind Sensitivity (Range) 
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Figure 3-22B Headwind Sensitivity (Range) 
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Figure 3-23 Tailwind Sensitivity (Range) 
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Figure 3-23B Tailwind Sensitivity (Range) 
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Figure 3-23C Tailwind Sensitivity (Range) 
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APPENDIX D PSEUDO-VACUUM TRAJECTORY PROGRAM 


DIM X(3ad),Y (302), RO(SAA), (Sag), MDOT (3aa),mM(3a2), Vi saa) 


DIM CDO (Sma) , MA (Siam) 
DIM T (Sia) 


BS5S9816S4000276F# 


LPRINT "TIME", "X","y", "MASS", "SPEED" 
FOR I=1 TO 32a 
X (1) =VO*COS (THETA) ¥T (1) 
¥ (1) =(VO*#SIN (THETA) *T (1) )— (16. WB7*T (1) “2D 
V (I) =( (VO*COS (THETA) ) “2+ (VO#SIN (THETA) — 
32. 174*T(1))*2)*.5 
TEMP =518. 688-. aaS5662*Y (1) 


Reto =e tien cee (i (ON s566eryY Cl) 7/Sts6. 668) )*S. 256) / 


(1716. S*TEMP) 
IF Y(I))=36089! THEN RO(I)=. 001271 
IF Y(1I))=36089! THEN TEMP=389. 99 
MA(I)=V(1)/SQR(1. 4*1716. 3* TEMP) 
CDO(1) =. 5@15-. 1925*MA (I) 
D(I) =. 276%. S¥RO(I) *. 3491*V(1) “2 
D(I)=D(1) /. 276*CDO(I) 
MDOT (I+1)=D(I)/ISP 
M(I+1)=M(1) -MDOT (I) *DELT 
T(I+1)=T(1)+DELT 
PSeMin oy Cl uaX (1), VOL), MiCl se vr) 
IF M(I+1) <=MMIN GOTO 22a 
NEXT I 


LPRINT "RAMJET BURNOUT BETWEEN T="5T(1I) 3 "AND"5T(I+1) 


Ww 
fia 


burning.Burnout will occur on the vacuum trajectory but not 


at the same point, as shown in Figure C-2. 


Figure C-2 Burnout Variations 


91 


In Figure C-1 wind (Vw) has been added for the curve #4 
in the x direction.The new equations are 


mx4 = (-Dx +) Dw) + Tx 


mZz4 


-Dz + Tz -mg 
where Dx + Dw= 1/2 9° Ug2a Ca A 


and Ux(wind relative to the projectile) =Vx+Vw. With no 
thrust curve 4 results. 
Finally, noe thrust. is,,added -to.requal , :drag. the 


following equations are valid: 


mMXS =O or Tx = Dy + Dw 
MZ5 = -mg or Tz = Dz 


The accelerometer(s) read the acceleration caused by any 
differences in the thrust and drag forces (zero in this 
case). 

ti east octear!, that! curves: InoandoSware equal Cand) 3» if 
thrust=drag) but the burntime of a ramjet will not be equal 


in each case due to different drag profiles. 


t = (mass fuel) (Isp)/thrust 


Therefore headwinds will require more thrust and 


burntime will be shorter.Conversely tailwinds allow longer 


174) 


case Ux=Vx and Uz=Vz.This trajectory is represented by curve 
2.An accelometer in the x direction measures Dx/m and Dz/mn 
in the z direction. 

If thrust is added, the trajectory is as shown at curve 


3(thrust greater than drag).The key equations are then: 


mX3 = -Dy ste Tx 
MZ g = =—D> + Te -mg, 
Accelerometers would measure (Tx-Dx)/m and (Ty-Dy)/m 


respectively.If thrust equals drag the trajectory would be 


curve 1. 


Figure C-1 Pseudo-Vacuum Trajectory Concept Curves 
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APPENDIX C-PSEUDO-VACUUM TRAJECTORY CONCEPT 


The basic concepts behind the pseudo-vacuum are not 
always readily apparent.This appendix analyzes the basic 
equations of motion in the following sequence of cases to 
illuatratea the concept: 

vacuum trajectory 
atmosphere trajectory 
thrusted trajectory 
wind influenced trajectory 
thrust equal drag trajectory 
In Figure C-1 (curve 1) the applicable equations of motion 


are: 


mX41 0 


ee 


M21 


-mg 


An accelometer on the projectile in the z direction would 
measure no acceleration as it is in free fall already. This 
is similar to a satellite orbit or someone on an elevator in 
free fall. 


When atmosphere is added the key equations are: 


| 
x 
N) 
1 


-Dy = -1/2 ° UL em aiCg 


=} 
N 
N 
" 


-Dz - mg 


where U is relative wind and D is atmospheric drag.In this 
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When t=0, the equation reduces to: 

-2Sw +1¢(28w) = 0 which is expected if the initial 
condition ‘ws zero. Gile. dispersion at zero equals zero at 
time zero). 

As t approaches infinity the equation value approaches 
infinite dispersion, but the oscillations of alpha have 
damped out. Note that the first derivative, proportional to 


velocity, tends to a constant value: 


Cw/ ($2 + w2 ))/m as t approaches infinity which is 


also to be expected. 
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APPENDIX B TRANSIENT ANGLE OF ATTACK CALCULATIONS 


The following equation of motion applies to the dynamics 
of angle of attack Calpha) oscillations due to the torque 


equation described in the three-degrees-of-freedom model. 


my = Lw which is proportional to alphact) 
where angle of attack is of the form 
-$t 
eo (t)=e sin(wt) where Sis the damping constant 


and w is the angular velocity. Therefore: 


+ 
. CCE dit) aT. 


(oy me) 


my 


a ne 
1 e-St sin(wt) dt dT 
079 


esol «-Ssintwl) i wlcoscwr).+ wl dT 


aie 


= wt(S§2 +w2) -28w + 


e-St(§2sin(wt) +2S8weos(wt) -w2sin(wt)) 


CiSeZ a ogee 
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impulse=950 “seconds” 

Mass Property Equations: 
CG=(225x39.4)-(52.O0xfuel burned)/(225-fuel burned) in 
IzzZ=.35-¢€.0019xSuel «burned? insslug-=tt 


Ixx=Iyy=18.8-(.0423xfuel burned) in slug-ft 


Aerodynamic Coefficient Equations: 
Cd=Cdo *.5 x Clex einGes) 
Cdo= =.1025 x MACH + 7.5025 (Los< MACH) < 2.2) 
Cm=Cl (Xcg-Xcp)/d 


Cl=Cl x a4 


Key Wind Tunnel Results (per degree): 
MACH=1.8 


Cdo=.317 Cmq=-203 


«x 6) zZ a 6 8 10 
cl @) -18 “38 sy oe das 
Xcp 43.34 45.6 48.2 “49.7 ("S022 Joo.0 
MACH=2.2 
Cdo=.276 Cmq=-190 
x '@) Pa 4 6 8 10 
Sa re) «LS .s40"%S3559"" 55 L238 


Xcp 42.9 42.9 44.6 45.1 46.0 44.6 
Cmq=d¢(Cm) /d(qd/2V) 


Initial Muzzle Spin Rate 30 PI radians /sec 


APPENDIX A-PROJECTILE CHARACTERISTICS 
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Figure A-1 ARMT Projectile Baseline Design 


The following characteristics are obtained from 


References 2 and 6. 

length=72 inches 

diameter=8 inches 

launch weight=225 lbm 

launch CG=39.4 inches from spike apex 


fuel weight=27 lbm (CG at 52.0 inches from apex) 


84 


except where harmonic content near resonance is 
contained in the wind profile. 

-headwinds create greater dispersion as ramjet 
burnout occurs before terminal guidance is 
initiated.Crosswinds and tailwinds generally will 


have minimal effects. 


Recommendations for further study related to this 
thesis are: 
-Dynamics. at . launch -inecluding-the-" initial spin 
damping, Magnus effects and aerodynamic jump. 
-Dynamics at terminal guidance including air inlet 
spike ejection and canard deployment. 
-Complete adaptation of BRL HTRAJ on the NPGS IBM 


system to use non-linear aerodynamics. 


isa) 
i) 


IV. CONCLUSIONS 

Generally this thesis confirma the viability of the 
pseudo-vacuum trajectory from. iah. Ehrust equals drag 
projectile.Indeed winds up to 120 feet per second and thrust 
errors up to five percent created dispersions of less than 
600 feet over a range of 175090 feet for the ARMT 
projectile. 

Specifically the following conclusions are made: 

-the drag due to wind (Dw) is the aerodynamic 
force which creates the most dispersion. 

-the net lift due to transients (Lw) is a minor 
aerodynamic force in creating dispersion. 

-the dispersion is proportional to the damping 
coefficient (Cmq),although for practical values of Cmq 
dispersion can be treated as insensitive to Cmq. 

-the projectile resonance or short period is at 
approximately 0). 22 Hz and is unstable in the 
resonance regime due to variable wind profiles.The phugoid 
is irrelevant. 

-thrust errors are the most significant errors 
contributing to dispersion. 

-winds have a minor effect on dispersion during 


powered fiigne® when thrust approximates drag, 


isa) 
fis 


Figure 3-28 Three Dimensional Trajectory With Altitude 
Dependent Wind Profile 
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Wind effects during near thrust equals drag portions of 
the trajectory were very small.The eritical factor 
affecting the dispersion will be the setting of the terminal 
guidance clock in the projectile.This presents several 
options for setting initiation of terminal guidance: 

-Fixed time irrespective of trajectory (ie. 
elevation of gun).Burn time would be shorter at less than 45 
degree launch and the basket would be more sensitive to wind 
and thrust errors.As well the trajectory would change and so 
would the predicted basket. 

-Fixed time for a given trajectory.This requires 
analysis of the sensitivity to elevation of the gun and a 
manual Tc setting method for gunners. 

-Ramjet burnout.This would ensure that Tc occurred 
at an “optimum"™ point almost exactly on the pseudo-vacuum 
Eaajectory, but the, Location along it at jiTe would be 
Saericult to predict. 

-At given pressure for altitude (Reference 5 ].This 
would establish accurate baskets for given gun elevations 
assuming the altitude chosen was before ramjet burnout. 

Several trajectories were simulated with a varying wind 
profile as) aa ftunction, of praltitude in the Standard 
Atmosphere.These are shown at Figure 3-28 as compared to the 


pseudo-vacuum trajectory. 
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Velocities at time 101 seconds were measured and are 
shown at Figure 3-27.The total variation of extreme cases is 


from 2285 feet/sec to 2315 feet/sec, or only 30 feet/sec. 
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Figure 3-27 Velocity Sensitivity 
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Deviations from the ideal trajectory were correlated 
and tabulated in Table 3.5.ResultS were conservatively 
linearized and are presented in that manner. Wind errors are 
meade in addition to the thrust errors at the left of the 
table.For example: 


If thrust error=-5% 
crosswind=120 ft/sec 


predicted dispersion is X=-540 ft 


Vo Loo it 
Z=130 ft ;from ideal basket. 


Table 3.5 Dispersion Analysis Results 


Thrust Errors Winds 
continuous 120 ft/sec 
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Figure 3-26 Tail/Headwind Effects (Burntime) 
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Figure 3-25C Tailwind Sensitivity (Height) 
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Figure 3-25B Tailwind Sensitivity (Height) 
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Figure 3-25A Tailwind Sensitivity (Height) 
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Figure 3-25 Tailwind’ Sensitivity 
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Figure 3-24C Headwind Sensitivity (Height) 
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Figure 3-24B Headwind Sensitivity (Height) 
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Figure 3-24A Headwind Sensitivity (Height) 
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Figure 3-24 Headwind Sensitivity (Height) 
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